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We report magnetoconductivity measurements on Y Ba^CuzO^ (x = 6.25 and 6.36) single crystals. Our main result 
is that both the in-plane Aa a b and out-of-plane Aa c magnetoconductivities exhibit the field dependence character- 
istic of "two-dimensional" quantum interference in applied magnetic fields H\\c. Namely, A<T Cjtl j, oc InH/Ho > 0, 
with Aa c /a c substantially greater than Aa a b/o- a b- We interpret this result as an evidence of interlayer incoherence 
in these crystals, so that the phase-coherent trajectories are mostly confined to one bilayer. 
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I. INTRODUCTION 

Understanding the anomalous features of the in- 
plane and out-of-plane normal state transport in lay- 
ered cuprates remains a challenge. Among the proper- 
ties of the resistivity tensor {p Cl p a h} that are some of the 
most unusual from the point of view of the conventional 
Fermi liquid theory are strongly temperature dependent 
anisotropy p c / p a b and coexistence of metallic p ab (T) and 
nonmetallic p c {T) over extended temperature and dop- 
ing ranges. It has been suggested that these anomalies of 
the out-of-plane transport are the result of interlayer in- 
coherence jjJH • This means that the phase coherent tra- 
jectories along which the single electrons maintain their 
phase memory are all confined to a single bilayer because 
the interbilayer transitions lead to dephasing of still un- 
known origin. An experimental confirmation of c-axis 
incoherence was provided by optical measurements 

The suggestion of interlayer incoherence in the nor- 
mal state of the cuprates, however, is not fully accepted. 
There are alternative models of the out-of-plane trans- 
port that do not involve interlayer incoherence QJ^] . The 
phase coherent trajectories in such models are three di- 
mensional and extend over many unit cells in both the 
in- and out-of-plane directions. 

Magnetoeffects due to quantum interference allow, in 
principle, the determination of the dimensionality of the 
phase coherence. Magnetoconductivity arises in rela- 
tively small magnetic fields due to contributions of the 
self-intersecting phase coherent trajectories along which 
the loops can be traversed in two different directions 
|^|,0J|] . The magnitude and field dependence of magne- 
toresistivity are determined by the probability of a tra- 
jectory to form a loop of a given area. If the phase- 
coherent trajectories are two-dimensional 2D, they are 
substantially more likely to form large loops than three- 
dimensional 3D ones. This results in a more pronounced 
effect, observable at higher temperatures. The field de- 



pendence of magnetoconductivity is also different in these 
two cases; namely, 2D trajectories lead to Act cx ln_ff, 
while 3D trajectories lead to Act cx H 1 / 2 . 

Here we present magnetoresistivities of two strongly 
underdoped Y BCL2CU3O5.25 and Y Ba2Cu^Oe.3e single 
crystals. Our results can be summarized as follows: Both 
components Ap c /p c and Ap a b/p a b of the magnetoresis- 
tivity tensor reveal the presence of two different conduc- 
tion mechanisms in magnetic fields H\\c. (1) In low fields, 
the magnetoresistivities are negative and have a logarith- 
mic field dependence over about a decade (at the low- 
est accessible temperature). This behavior represent the 
main focus of this paper. What distinguishes these re- 
sults from traditional experiments on ultrathin films and 
other 2D systems ||[)), is that the out-of-plane magne- 
toresistivity has the logarithmic field dependence and is 
substantially larger than the in-plane magnetoresistiv- 
ity. We interpret these findings as a confirmation of a 
strong interbilayer decoherence and, therefore, the two- 
dimensional nature of the phase-coherent paths. Under 
this condition, both components of the conductivity de- 
pend on the in-plane phase coherence length £ v , which 
changes with field logarithmically, thus producing the 
corresponding independence of the magnetoresistivities. 
(2) In high magnetic fields and at low temperatures, both 
magnetoresistivities become positive and change as jH 2 
up to the largest available field of 14 T. This contribution 
arises from antiferromagnetic correlations. 

II. EXPERIMENTAL DETAILS 

Strongly underdoped Y Ba^Cu-^Ox (x = 6.36 and 6.25) 
single crystals were grown in gold crucibles using the self- 
flux method. The oxygen stoichiometry was adjusted 
by annealing the samples at 500°C in a predetermined 
O2—N2 atmosphere [ ^0[ . In-plane p a b and out-of-plane p c 
resistivities and the respective magnetoresistivities (MR) 
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were measured by a multiterminal method on the same 
single crystal, in magnetic fields H up to 14 T. This al- 
lows us to carry out a quantitative comparison between 
in-plane Ap a b/p a b and out-of-plane Ap c / p c magnetore- 
sistivities as a function of temperature T and applied 
magnetic field H . All previous measurements of MR in 
Y Ba2Cu^O x reported in the literature were obtained by 
a four-point method, so that Ap a b/p a b and Ap c / p c were 
measured on different single crystals. Since magnetore- 
sistivity is very small, inevitable variations of stoichiom- 
etry, shape, and size of specimens have a large effect on 
its value and preclude any quantitative analysis of corre- 
lations between different components of the magnetore- 
sistivity tensor. 

Low contact resistance gold leads were attached us- 
ing thermally treated silver pads and room tempera- 
ture silver epoxy. We applied an ac electrical current 
I ~ 0.1 mA through the two leads on the "top" face 
and measured the two voltage drops on the same (top 
voltage Vt) and the opposite (bottom voltage V b ) faces of 
the sample [jll). The two voltages Vt,& and their varia- 
tion with field SVt.b were measured using a low-frequency 
17 Hz resistance bridge at constant temperature, while 
sweeping the magnetic field H applied parallel to the c- 
axis of the single crystal. 

The temperature of the sample chamber was kept con- 
stant using two temperature sensors (Pt thermometer for 
T > 100 K and Cx sensor for T < 100 K). The tempera- 
ture change 5T{H) resulting from the magnetoresistance 
of the sensors was measured by carefully calibrating them 
using a capacitance thermometer, which has zero magne- 
toresistance. The effect of the magnetic field on Vt, b was 
then determined by subtracting the effect of this temper- 
ature change from the raw data; i.e., 

AV t , b (H) = 8V ttb -8T(H)^±. (1) 

Unlike in the four-point method, the distribution of 
the current inside the sample is nonuniform. As a result, 
the Hall voltage contributes to the field variation of the 
top and bottom voltages. To first order in magnetocon- 
ductivities, appropriate for these measurements, 

AV t ,b{H) = A t , b Aa xx + B tjb Aa zz + C t 

,b&xy i (2) 

where Aa xx and Aa zz are magnetoconductivities, while 
a xy is the off-diagonal component of the conductivity ten- 
sor. The coefficients A t . b , B t ^ and Ct,b can be obtained 
from the solution of Laplace's equation determining the 
distribution of voltage inside the sample. (Within a linear 
approximation, these coefficients depend on the zero-field 
values of the conductivity tensor, which is diagonal.) 

The most reliable way to determine the magnetocon- 
ductivities Aa xx and Aa zz from Eq. (2) is to elim- 
inate the contribution of the Hall effect by measur- 
ing AVt j b(iJ) at two opposite polarities ±H of the ap- 
plied magnetic field. The even component AV t + b (H) = 



l/2[AV t , b (H) + AV t ,b{-H)] determines the magneto- 
conductivity, while the odd component AV t ~ b (H) = 
l/2[AV t , b (H) - AVt, b (-H)] = C t .bO xy . The magneto- 
conductivities were then calculated from AV t + b (H) using 
an algorithm described in O] . 

As an example of the raw data, the main panel of Fig. 
1 shows the field dependence of AVt of the x — 6.36 
single crystal, measured at T = 100 K, while the inset 
displays its odd AV t ~ and even AV t + components. The 
field dependence of AVt clearly separates into quadratic 
dependence of AV t + and, consequently, of magnetocon- 
ductivity, and linear dependence of AV t ~ and, hence, Hall 
conductivity a xy . 



20 




H(T) 

FIG. 1. Main panel: The field H dependence of 
the top voltage AV t for the x = 6.36 single crys- 
tal, measured at T = 100 K in H c; In- 
set: odd AV t ~ = l/2[AV t (H) - AV t (-H)} and even 
AV t + = l/2[AVt(H) + AV t (-H)] components of the top volt- 
age obtained from the data shown in the main panel. 

The multiterminal technique has never been used be- 
fore to measure magnetoconductivities of cuprates. The 
application of this method to magneto-measurements has 
subtle, but important, differences with the conventional 
four-point method. For example, in the multiterminal 
method one measures magnetoconductivities, while in 
the conventional four-point method one measures magne- 
toresistivities. It should be noted that most theoretical 
approaches to magnetoeffects describe magnetoconduc- 
tivities rather than magnetoresistivities. Assuming that 
a xz and a yz are negligible, the relationship between mag- 
netoresistivities and magnetoconductivities is given by 
& 

^Pxx ^&xx ®xy 

Ap zz Aa zz 

Pxx @xx ® xx Pzz O 'zz 

Comparing the values of the even and odd components 
AV^, we find that in most cases Aa xx ~ a xy and, 
since both of them are very small, the quadratic term 
in Eq. (3) is negligible. Therefore, the magnetoresis- 
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tivities are given directly by magnetoconductivities with 
inverted sign. The only exception to this case is at very 
low magnetic fields where a xy > Aa xx since a xy oc H 
while Aa xx oc H 2 . Therefore, at very low fields, one 
needs to know a xy in order to obtain magnetoresistivity 
from magnetoconductivity. In this paper, we present all 
our results in terms of magnetoresistivity obtained from 
the measured magnetoconductivity according to Eq. (3) 
where we neglect the quadratic term. This facilitates 
a direct comparison between our data and already pub- 
lished magnetoresistivity data. 



III. MAGNETORESISTIVITY 

Figures 2(a) and 2(b) show the T-dependences of the 
resistivities of Y Ba2Cu3Oe.se and Y Ba2CusOe.25 sin- 
gle crystals, respectively, in zero magnetic field. These 
single crystals display a diverse T— dependence of resis- 
tivities, and antifcrromagnctism below T i=s 40 K in 
x = 6.36 sample and for all T < 300 K in x = 6.25 
crystal. Coexistence of metallic p a b and nonmetallic p c 
over an extensive temperature range is characteristic of 
underdoped cuprates. We note that the in-plane resis- 
tivity p a b is metallic above a certain temperature even 
in the strongly underdoped samples. The crossover to 
nonmetallic behavior in p a b takes place at progressively 
higher temperatures in more underdoped samples. For 
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FIG. 2. Temperature T dependence of zero-field in-plane 
Pab and out-of-plane p c resistivities for (a) Y Ba2CuzOe.M>, 
and (b) Y BaiCu-sOti.2h- 



example, the minimum in p a b shifts from w 50 K for the 
x = 6.36 sample to ~ 175 K for the x = 6.25 sample. 
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FIG. 3. Field H dependence of magnetoresistivities for 
H || c: (a) Apab/pab and (b) Ap c / p c for YBa2Cu 3 06.36 at 
50 K and 75 K. 

Measurements of the MR tensor on YBa2Cu3Oe.se at 
temperatures higher than 75 K showed that the field 
dependence of both in-plane and out-of-plane MR is 
quadratic However, at lower temperatures, the field 
dependence of the MR tensor becomes non trivial [see 
Figs. 3(a) and 3(b)]. Namely, at 50 and 75 K, a nega- 
tive MR component, with a large slope at relatively low 
fields, is superimposed on a positive quadratic field de- 
pendence. Notice that the minimum in Ap c /p c is ap- 
proximately twice the minimum in Ap a b/p a b- At higher 
magnetic fields, MR is positive and changes quadratically 
with H . 

Figures 4(a) and 4(b) display the magnetic field H 
dependence of MR of the Y Ba2Cu^Oe.ih single crystal. 
The same phenomenon as in the previous sample, namely, 
sharply developing negative MR in low fields which turns 
into positive quadratic dependence in high fields, is even 
more pronounced. 

The field dependence of MR illustrated in Figs. 3 and 
4 can be described for both samples as a sum of two 
contributions: a negative contribution Qi{H) where i = 
{c, ab}, and a positive quadratic contribution; i.e., 



3 



^ = Qi(H) 

Pi 



jiH 2 ; Qi(H) <0; 7l >0. 



(4) 



We examine in more detail the field dependence of the 
MR components measured at 75 K on Y Bcl-zCuzOq^ 
in Fig. 5 (a). First, we note that the negative contri- 
bution to MR is much larger for the out-of-plane com- 
ponent than for the in-plane component; for example, 
the absolute value of Ap c / p c at its minimum is about 7 
times the corresponding Ap a b/p a b- Second, both Ap c /p c 
and Ap a b/pab measured in high magnetic fields are well 
fitted with a parabolic dependence jiH 2 — e;, as shown 
by the solid curves, with j c = 1.3 x 10~ 5 T~ 2 , j a b — 
6.25 x 1(T 6 T- 2 , e c » 0.15%, and e ab » 0.03%. We dis- 
cuss the origin of this quadratic H- dependence at the 
end of the next Section. 
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FIG. 4. Field H dependence of magnetoresistivities: (a) 
Ap^/pab and (b) Ap c /p c for YBa 2 Cu 3 6 .25- 

The field dependence of Qi(H) in Eq. (4) can be de- 
termined by subtracting ^iH 2 from the total MR, and is 
illustrated in the semi-log plot of Fig. 5(b). The straight 
lines are fits to the data and reveal the logarithmic de- 
pendence of Qi(H) within a certain range of H. The 
functional dependence of Qi{H) can be summarized as 
follows: 



-€i H 2 /H 2 

ln(g/g Q ) 
~ tl ln(Hi/ff ) 



H<H 
H < H < Hi 
H>H X 



(5) 



covering almost a decade of H . The value of the parame- 
ter Hq is determined by the intercept of the straight line 
according to Eq. (5). For H > H x « 4 T, Qi(H) satu- 
rates at the value of obtained earlier by the parabolic 
fit in Fig. 5(a). The \n(H\/ Hq) in the denominator of 
Eq. (5) is required to match the interpolations of Qi(H) 
at H < H < Hi and H > Hi. 
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The most interesting part is the logarithmic dependence, 



FIG. 5. (a) Magnetoresistivities Ap c /p c and Ap a b/p a b for 
YBa 2 Cu 3 6 .25 at T = 75 K in magnetic fields H\\c . The 
solid lines represent quadratic fits ■jiH 2 — €i (i = {c, ab}). (b) 
Orbital component of magnetoresistivity Qi — Api/pi —jiH 2 
plotted vs logarithm of field. The straight lines are guides to 
the eye corresponding to the logarithmic dependence given by 
Eq. (5). The crossover fields Ho and Hi are indicated. 

The negative contribution to the in-plane MR is 
smaller and, therefore, the data are noisier, but a sim- 
ilar procedure of fitting the high field regime data with a 
parabolic dependence yields Q a b(H) with approximately 
the same H-dependence as Q C (H); i.e., Q c {H)/Q ab (H) « 
const = 7. 



IV. DISCUSSION 

The most intriguing observation that emerges from the 
data presented above (Figs. 3 and 4) is the presence of a 
negative MR component in low fields [given by Qi (H) in 
Eq. (4)] , especially pronounced in Ap c /p c , which exhibits 
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at the lowest accessible temperatures a logarithmic in- 
dependence. This negative MR is superimposed on a 
positive quadratic background ["fiH 2 in Eq. (4)], which 
is due to antiferromagnetic ordering Jl4| . Below, we show 
that the low field features of MR are direct consequences 
of interlayer incoherence. 

Our approach is based on the understanding that 
strongly underdoped layered crystals like YBa 2 Cu^O x 
are characterized by incoherent out-of- plane transport. 
This was shown in optical measurements of the conduc- 
tivity on underdoped Y Ba 2 Cu-sO x J^]. In incoherent 
crystals, the value of the out-of-plane coherence length 
i VlC is T- independent, equal to its minimum possible 
value. Since the coherence length of electron wave func- 
tions in a crystal cannot be smaller than the size of the 
atomic orbitals that overlap over the interbilayer dis- 
tance, the minimum possible value of £ VyC is the inter- 
bilayer distance; i.e., £ ViC — £q — 11. 7A. Thus, a fun- 
damental property of fully incoherent crystals is that 
£ VtC does not change with temperature or magnetic field. 
Therefore, the only length scale which determines the 
dissipation and can change with temperature or applied 
magnetic field is the in-plane phase coherence length £ tp . 
(Hereafter, we assume for brevity that the planes are 
isotropic and omit the subscript ab.) Under these con- 
ditions, both conductivities depend only on the variable 
£<p [&ab{£tp) and <7 c {£tp)], so that their temperature and 
field dependences come from that of £ tp . (An instructive 
discussion regarding the dependence of conductivity on 
the spatial length scale of inelastic scattering processes 
is presented in ||). 

The immediate consequence is a correlation between 
the magnetoconductivities and the field variation of the 
in-plane coherence length at constant T: 

Aa ab (H) M V (H) Aa c (H) A£ V (H) 

a , = Kab — 5 — = Kc ~£ ■ (6) 

where re, = d In cr^ jd In £ v (i = {ab, c}). 

The coefficients dimensionless numbers of the 

order of unity. Their sign determines the type of electri- 
cal conduction (metallic or nonmetallic); i.e., if «j > 0, 
then daijdT < and the conduction is metallic, while if 
Ki < 0, then d<Ji/dT > and the conduction is nonmetal- 
lic Jl5[ |. It is important to remember that in Eq. (6) and 
hereafter, Aa a b/^ab and Aa c /a c are the "orbital" con- 
tributions to magnetoconductivity also denoted as Qi in 
Eq. (4). 

According to Eq. (6), the sign of the MR tensor is 
determined by the sign of k, and by the effect of the 
magnetic field on £ ip . It is known ||, and in the Ap- 
pendix shown in more detail, that the magnetic field 
decreases the phase coherence length due to its effect 
on self-intersecting trajectories, so that A£ ip (H) < 
in a weak field. The sign of Ki can be inferred from 
the T-dependence of the resistivity. Specifically, for 



YBa 2 Cu 3 6 .25, K c , ab {T) < for T < 175 K since both 
p c and p a b are nonmetallic [see Figs. 2(a) and 2(b)]. 
Then, according to Eq. (6), both magnetoconductivity 
components are positive (MR, correspondingly, are neg- 
ative) for T < 175 K, as indeed is the case [see Fig. 4]. 
By applying the same reasoning to the YBa2Cu 3 OQ.z& 
sample, one concludes that Ap c /p c should be negative, 
while Ap a b/pab should be practically zero at 75 K and 
should become negative for T < 50 K , where p a b turns 
nonmetallic as well. This is consistent with the low field 
data of Fig. 3. 

Also, according to Eq. (6), the field dependence of 
both MR is the same, determined by that of A£ V (H). 
Therefore, their ratio should be a constant, given by the 
ratio of the corresponding k%; i.e., 

A(ln cr c ) _ Kc_ 
A(lner afc ) K ab ' 

As discussed in Refs. |il^ , |l7j ], in incoherent crystals 
a c (£ v ) oc <J a b{£^)/£%- As a result, k c = K ab - 2. There- 
fore, when K a b is negative (as is the case with our crystals 
at low temperatures), the absolute value of A (In a c ) is 
greater than that of A(ln a a b), namely: 

AjWcj = Kb\ +2 ^ 
A(ln a a b) | Kab | 

In other words, when the T-dependence of both conduc- 
tivities is nonmetallic, the out-of-plane conductivity de- 
pends stronger on the phase coherence length, and the 
respective MR is greater by a constant factor. 

This conclusion is, indeed, supported by the exper- 
imental data. The ratio of the magnetoresistivities is 
a constant factor which can be estimated, for example, 
from Fig. 5(b). Namely, for YBa 2 Cu 3 6 .25 at T = 75 K, 
Q c is about seven times greater than Q a b- From Eq. (8), 
we estimate n a b ~ —0.33. 

We turn now our attention to the specific field depen- 
dence of the magnetoresistivities. The logarithmic field 
dependence of Q at low temperatures and in relatively 
small magnetic fields indicates 2D quantum interference 
HH . This effect is due to the presence of self- intersecting 
trajectories along which electrons can traverse the loop 
in opposite directions (6|^]. A small applied magnetic 
field gradually destroys the quantum interference, which 
induces a variation in the phase coherence length and, 
therefore, a variation in the conductivity. 

The self-intersecting trajectories constitute a small 
fraction of the total number of phase coherent trajec- 
tories of a given length. The great majority of the phase 
coherent trajectories (trajectories over which the phase 
changes predictably and, therefore, no irreversible pro- 
cesses are involved) are not self-intersecting. 

As Fig. 2 shows, both resistivities change strongly 
with temperature and, therefore, their magnitude and 
T-dependence are determined by the contribution of the 
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majority, non-intersecting trajectories. The small, tem- 
perature dependent corrections to conductivities (in zero 
field) due to self-intersecting trajectories, known as cor- 
rections due to weak localization, can only be observed 
when the main contribution to conductivity which is due 
to non-intersecting paths saturates (regime of residual 
resistivity). Our samples are not in this regime. There- 
fore, the contribution of the trajectories with loops to the 
T-dependence of both conductivities is negligible in our 
samples and a a b(T) and <J C {T) do not reflect the charac- 
teristic InT (in 2D) or T 1 / 2 - T 1 / 3 (in 3D) dependence 
due to quantum interference, which is usually the subject 
of discussion in the literature. 

On the other hand, the magnetoconductivity produced 
by a weak magnetic field is due to self-intersecting trajec- 
tories because all the other (conventional) contributions 
to magnetoconductivity are still negligible. Therefore, 
the magnetoconductivity caused by the effect of the field 
on the self-intersecting trajectories can be experimentally 
identified. 

In order to determine the field dependence of both 
magnetoconductivities, we need to evaluate the magnetic 
field variation of the coherence length \A£ V (H) / £ ip \. A 
semi-quantitative derivation of this dependence is given 
in the Appendix, with the following results: 
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-V H 2 /H 2 

\n(H/H„) 



H<H 
H <H <H X 
H > H x 



(9) 



Here rj is a small number determined by the relative 
weight of the number of the self-intersecting trajectories 
with respect to that of the majority nonintersecting tra- 
jectories. The lower crossover field Hq is determined by 
the condition that the field flux through the largest loops 
along the phase coherent trajectories is approximately 
equal to the flux quantum <fio = 2 x 10 _7 Oe cm 2 ; namely 
Ho ~ 0o It- %■ The upper crossover field Hi is determined 
by the similar condition that the field flux is of the order 
of <fio through the smallest possible loops (of the order 
of the elastic mean free path l e i); i.e., Hi ~ 4>q/1 2 i- For 
H > Hi , no trajectory with loops contributes to the value 
of the phase coherence length. Hence, A£ v saturates. 

From Eqs. (6) and (9) follows Eq. (5) with e$ = —rj Ki. 
The nontrivial logarithmic field dependence appears only 
when there is a significant field range between Hq and 
Hi . This requires the phase coherence length to be sub- 
stantially greater than the elastic mean free path, which 
means a well developed diffusive regime. Therefore, the 
magnetoeffects due to self-intersecting trajectories are 
more pronounced in stronger underdoped samples. For 
YBa 2 Cu 3 6 .2 5 at T = 75 K [see Fig. 5(b)], we esti- 
mate that e c « 0.15%, Hq ~ 0.3 T which corresponds to 
l v ~ 10 3 A, and Hi ~ 4 T (£ el ~ 270 A). For H > 4 T, 
the contribution of quantum interference to magnetocon- 
ductivities saturates and other types of contributions to 



MR become dominant, as discussed below. 

Experimentally, one can relatively easy distinguish 
between the cases of two-dimensional and three- 
dimensional phase coherent trajectories. First, the mag- 
netoeffect of self-intersecting trajectories is substantially 
greater in 2D, because of the greater probability of a 2D 
trajectory to form a large loop (77 ~ \/£ e i in 2D and 
77 ~ A 2 /£ 2 el in 3D, see the Appendix). Therefore, in this 
case, the effect becomes pronounced and observable at 
relatively high temperatures. Second, the logarithmic H 
dependence can be distinguished from H 1 ! 2 rather clearly 
even for only one decade in H . These features observed 
in MR would point to the 2D nature of the phase co- 
herent trajectories. The two manifested together, as is 
the case with our MR data of Y Ba2Cu%O x shown above, 
make this conclusion very compelling. It should be noted 
that the "three dimensional" H 1 / 2 dependence of MR in 
manganites was reported recently by QingAn Li, Gray, 
and Mitchell g||. 

Finally, we discuss briefly the origin of the positive con- 
tribution that changes as jiH 2 and becomes dominant at 
applied magnetic fields larger than 4 T (see Figs. 3 and 
4) . A similar type of MR was reported by other groups 
as well 14 1 , and it was associated with the antiferromag- 
netic fluctuations and antiferromagnetic ordering. In- 
deed, YBo,2Cu30q,25 is antiferromagnetic for all temper- 
atures below room temperature, while Y Ba2CusOe.36 is 
antiferromagnetic for T < 40 K . In high magnetic fields, 
the contribution of the AF correlations to MR dominates 
the smaller contribution of the self-intersecting trajecto- 
ries. At sufficiently low temperatures, these two contri- 
butions can be clearly separated as illustrated in Figs. 
5(a) and 5(b). 



V. SUMMARY 

We present magnetoresistivity data for two strongly 
underdoped single crystals oiYBa2Cu^O x with x = 6.36 
and 6.25. Both in-plane Ap a b and out-of-plane Ap c mag- 
netoresistivities MR were measured simultaneously on 
the same single crystal using a multiterminal method. 
This permits a direct comparison between their temper- 
ature and field dependences. 

The most interesting observation is a negative MR in 
low applied magnetic fields. This low field contribution is 
characterized by two important features: First, the effect 
is strong and, therefore, is pronounced even at relatively 
high temperatures (75 K and higher). Second, the field 
dependence is consistent with \n(H), rather than H 1 ! 2 . 
These features point towards the two dimensional na- 
ture of the phase coherence in these crystals; namely, the 
phase-coherent volume contains only one or two neigh- 
boring bilayers, while the coherence length along the 
planes is orders of magnitude greater than the size of 
the unit cell. 
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We attribute the second contribution to MR, which 
changes quadratically with field up to H = 14 T, to 
the effect of H on the AF correlations. This positive 
contribution dominates at fields above 4 T and induces a 
change in the sign of the total MR for both components 
of the magnetoresistivity. 



VI. APPENDIX 

To determine the effect of the magnetic field on self- 
intersecting phase coherent trajectories, and the resul- 
tant effect on the phase coherence length, we use the 
qualitative approach described in §3,0]. The coherent 
trajectories consist of two groups. The great majority 
of them are not self-intersecting and they determine the 
value and the temperature dependence of the phase co- 
herence length and, consequently, of the conductivities. 
A much smaller fraction of the coherent trajectories have 
loops and exhibit the interference effect. 

We begin our estimate by defining the phase coherence 
length in zero magnetic field £ v (0) as the mean square av- 
erage distance in the x-direction that an electron travels 
while retaining its phase memory (i.e., without encoun- 
tering dephasing inelastic collisions): 



4(0) 



x 2 [(l~r 1 )P Q (x)+r 1 P sl (x)}dx. 



(10) 



Here Po(x) [P S i(x)} is the probability that the phase co- 
herence is retained between points A and B separated 
by a distance x, when all trajectories without [with] self- 
intersections are counted. Both probabilities are nor- 
malized to unity. The small number 7/ <C 1 reflects the 
relative weight of self-intersecting trajectories. We also 
introduce p(x, A) as the probability density that a self- 
intersecting, coherent trajectory (in zero field) has a loop 
of area A, subject to the normalization condition: 



/ p[x,A)dA = P si {x). 
UP, 



(11) 



Here, the smallest possible area of the loop is determined 
by the elastic mean free path £ e i , and the largest possible 
loop is determined by the distance x between A and B. 

Along self-intersecting trajectories, an electron can tra- 
verse the loop in opposite directions with no phase dif- 
ference between the respective amplitudes. An applied 
magnetic field, however, introduces a phase difference 
ipH = 27r0/0o between the two alternative routes along 
the same trajectory. Here 4> — AH is the magnetic field 
flux through the loop, A is the area of the loop normal to 
the field, and 4>q = 2 X 10~ 7 Oe cm 2 is the flux quantum. 
As a result, when an electron arrives at point B, follow- 
ing a self-intersecting route, its phase is unpredictable, 
and varies by ifn > 7r if the corresponding trajectory has 



a loop with area A > 4>q/2H. Therefore, such a trajec- 
tory is no longer phase coherent and does not contribute 
to the average given by Eq. (10). This is the mechanism 
by which a weak field changes the coherence length and, 
respectively, the conductivities. 

We give a rough estimate of A£ V (H) by considering 
that all sclf-intcrsecting trajectories with loop areas A > 
4>o/2H do not contribute to the average in Eq. (10), while 
those with A < <fio/2H still do. The phase coherence 
length in a magnetic field is then given by: 



x dx 



{l-r))P (x) + ri 



p(x,A)dA 



(12) 

where the magnetic length £fj — (Hc/2eH) 1 / 2 is deter- 
mined by the condition ir£ 2 H = (f>o/2H. 

Taking into account Eqs. (10) and (11) and the small 
value of ry, we obtain from Eq. (12): 



2A£ V (H) _ J^ x 2 dxJ^ H p(x,A)dA 
J °° x 2 Po(x)dx 



-T) 



(13) 



Thus, the decrease of the phase coherence length A£ V (H) 
is proportional to the weight of excluded self-intersecting 
trajectories, namely, those with loop areas A > Tt£ 2 H . 
The lower limit of integration in x is £h because the 
trajectories with a coherent distance x < £h cannot have 
loops with area A > ir£ 2 H and, therefore, are not part of 
the excluded self-intersecting trajectories. 

Next we find the expression for p(x, A) using the dif- 
fusive approximation H|. Since phase-coherent trajec- 
tories are two-dimensional (we discuss here incoherent 
crystals), the probability to find an electron at a point 
r from its starting origin at t = 0, is W(r)d r cx 
cxp(—r 2 /4Dt)d 2 r/4Dt, where D is the difussion coeffi- 
cient. The trajectory intersects itself between the times 
t and t + dt if it enters the area Xvdt around the origin. 
Here A is the de Broigle wavelength (the trajectory is 
viewed as a tape of width A, rather than a line). The 
probability of return, and, consequently, the probability 
to form a loop is given by p(t)dt cx W(0)Xvdt oc Xvdt/Dt. 
Since the area of the loop is proportional to the average 
distance from the origin, A oc (r 2 ) oc Dt, and dA oc Ddt, 
the probability of a loop with area between A and A+dA 
is given by (Xv/D)dA/A. Thus, the probability that a 
2D trajectory has a loop of area A scales as 1/A. 

The prefactor Xv/D determines the overall probability 
of self-intersecting trajectories and gives the number r), 
which wc introduced "by hand" earlier: 77 ~ Xv/D ~ 
X/£ e i, because D ~ v£ e i. Assuming A ~ 2 — 3 A 
and £ e i ~ 200 — 300 A, which is a reasonable estimate 
given the density of charge carriers, we obtain r\ ~ 0.01. 
ft means that the maximum effect of removing self- 
intersecting trajectories from the average in Eq. (13) re- 
duces £ v by about one percent. Other numerical factors 
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also absorbed in rj can reduce it further. This correlates 
with the maximum value of the negative MR in Fig. 5(b) 
of the order of 0.1%. 

Since the diffusive motion has no memory, having a 
loop does not affect the statistical properties of the rest 
of the trajectory. Therefore, the probability density 
p(x, A)dA can be roughly estimated as: 

P sl {x) dA 



p{x,A)dA 



In(xV^) A 



(14) 



where the logarithm in the denominator is the normal- 
ization factor due to Eq. (11). Now the integral in (13) 
becomes: 



x dx 



p{x,A)dA ■ 



x 2 P si {x) - ) J.St dx. 



ln(*7*d) 



(15) 



The characteristic logarithmic field dependence of this 
integral appears when £ 2 H decreases below the value of 
l 2 v . This is equivalent to the condition that the field 
flux through the largest loops along the phase coherent 
trajectories is approximately equal to (f>o/2; hence, the 
characteristic field Hq = (/)q/2tt£ 2 ^. As £ v increases with 
decreasing temperature, the onset field Hq can become 
very small and the MR due to the interference effect man- 
ifests itself when all the other sources of MR are negligi- 
ble. The upper crossover field is defined by the condition 
that the field flux is 4>o/2 through the smallest loops with 
the area of the order of tt£ 2 i', hence Hi = (j>o/2nf^ l . The 
upper field does not change with temperature, while the 
lower field Hq decreases with decreasing T. Therefore, 
this effect is observable at sufficiently low temperatures 
when H < H < Hi or, equivalently, £ 2 > l\ > £\ v 

In this regime, the integral (15) can be easily estimated 
because the logarithms are slowly changing functions and 
they can be taken at the value x ~ £ v , which corresponds 
to the maximum in x 2 P S i(x). Other contributions to this 
integral are negligible in comparison with the logarith- 
mic increase ~ ln(^/^) = \d.{H/Hq). Correspondingly, 
the logarithm in the denominator of Eq. (15) becomes 
ln(£ 2 /£^i) = ln(Hi / Hq) . In this regime, the denominator 
in Eq. (13), J °° x 2 P Q (x)dx ~ J °° x 2 P si (x)dx « £ 2 ip , and 
we get Eq. (9), with all numerical factors absorbed in rj. 
Hence, according to Eq. (6), both conductivities have the 
logarithmic field dependence. For H > Hi, the integral 
(15), hence, Eq. (13) saturates since no self-intersecting 
trajectories contribute to the average (12). 

It is important to underline that the transport in these 
crystals is obviously three-dimensional, but the phase co- 
herence is two-dimensional, meaning that the distribu- 
tion of the area of the loops along phase coherent tra- 
jectories scales as 1/A, leading to the logarithmic field 
dependence of A£ lp (H) and both magnetoconductivities. 
This result remains valid even when the coherence is es- 
tablished between two neighboring bilayers. As long as 



£p,c ~ *C £tp,abi the distribution of loop areas is close 
to 1/A. 

Truly 3D phase coherent trajectories requires £ VtC 3> 
£q. Then, a similar analysis 



shows that the distribu- 
in the regime Hq <C H -C Hi, Eq. (A) gives 



tion of the loop areas is dA/A^I 2 and i] ~ \ 2 /£ 2 el - Then, 



A£ V {H) 



e 2 



H 

Hi 



1/2 



(16) 
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